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SUBSEA  PERMAFROST 
IN  HARRISON  BAY,  ALASKA 
An  Interpretation  from  Seismic  Data 


K.  Gerard  Neave  and  Paul  V.  Sellmann 


INTRODUCTION 

The  objective  of  this  study  was  to  obtain  infor¬ 
mation,  through  analysis  of  seismic  records,  on 
the  distribution  and  properties  of  subsea  perma¬ 
frost  in  the  Harrison  Bay  region  of  the  Alaskan 
Beaufort  Sea. 

Information  on  this  region  is  based  primarily  on 
seismic  studies.  No  core  drilling  or  sample  analysis 
has  been  done,  although  temperature  data  and  ob¬ 
servations  of  penetration  resistance  have  been  ob¬ 
tained  from  shallow  probe  observations  (Harrison 
and  Osterkamp  1981).  Several  high-resolution 
seismic  lines  were  run  by  Rogers  and  Morack 
(1981).  However,  the  penetration  achieved  during 
their  study  (40  m)  largely  precluded  the  detection 
of  deeper  high-velocity  material.  The  lack  of 
deeper  drilling  and  core  analysis  has  resulted  in  a 
critical  deficiency  in  ground  truth  data  to  verify 
the  seismic  interpretation. 

This  investigation  is  based  on  the  interpretation 
of  monitor  seismic  records  from  petroleum  ex¬ 
ploration  programs.  The  first-return  data  used  in¬ 
cluded  speculative  and  nonproprietary  data  from 
Western  Geophysical  Company  of  America,  plus 
data  released  by  British  Petroleum.  A  total  of  ap¬ 
proximately  450  km  of  seismic  line,  with  one 
record  per  kilometer,  was  examined.  The  locations 
of  these  lines  are  shown  in  Figure  1 . 

The  data  for  shoreline  transitions  and  other 
shallow  water  coastal  areas  are  from  a  winter  sur¬ 
vey  conducted  over  the  ice.  The  remaining  shot 
lines  are  from  a  marine  survey.  A  summary  of  the 
field  parameters  used  in  collecting  both  sets  of 
data  is  presented  in  Table  1. 

In  seismic  data  processing  for  petroleum  explor¬ 


ation  the  emphasis  is  normally  on  deep  targets, 
commonly  with  first-break  suppression,  increased 
gain  with  time,  time-variable  filters,  and  normal 
move-out  corrections.  These  procedures  tend  to 
compromise  the  quality  of  the  data  from  near  the 
surface.  To  obtain  as  much  information  as  possi¬ 
ble  from  the  records  without  costly  processing, 
monitor  records  were  produced  by  playing  back 
approximately  the  first  2  seconds  of  the  field  tapes 
with  expanded  gain.  They  were  printed  in  a  “wig- 
gletrace”  or  variable  area  format,  without  normal 
move-out  corrections,  or  filtering. 

An  appreciable  increase  in  seismic  velocity  oc¬ 
curs  when  most  unconsolidated  materials  freeze 
(Aptikaev  1964).  This  increase,  and  the  existence 
of  large  amounts  of  seismic  data  from  surveys 
conducted  for  petroleum  exploration,  make  seis¬ 
mic  techniques  a  reasonable  approach  for  investi¬ 
gating  the  distribution  of  ice-bonded  subsea  per¬ 
mafrost.  Hunter  et  al.  (1976)  and  Sellmann  et  al. 
(1980)  have  shown  that  when  records  are  avail¬ 
able,  and  their  quality  and  field  recording  param¬ 
eters  are  appropriate,  permafrost  data  can  be  ex¬ 
tracted  from  the  monitor  records.  Direct  wave  vel¬ 
ocities  and  refraction  interpretation  of  the  first  re¬ 
turns  can  give  information  on  velocity  structure 
(Fig.  2).  These  velocity  data  are  used  to  predict  the 
distribution  of  ice-bonded  permafrost  and,  when 
resolution  is  adequate,  the  depth  to  the  top  of  the 
frozen  sediments.  Reflection  interpretations  also 
supplement  the  depth  information  when  shallow 
reflectors  are  strong  enough. 

A  number  of  error  estimates  are  presented  in  the 
Spatial  Resolution  section  to  verify  that  the  meas¬ 
urements  and  analytical  methods  are  reasonable. 
These  estimates  are  not  intended  to  set  strict  limits 


Figure  1.  Harrison  Bay  area.  Solid  lines  are  from  a  marine  survey  and  dashed  tines  from  an  ice  survey. 


Table  1.  Harrison  Bay  survey  specifications. 


BP  Alaska  ice  survey 

Western  Geophysical  shallow  water  marine 

Source 

Explosives  (23  kg  or  SO  lb) 

10  air  guns,  9500  cm'  (580  in’.),  each  at  a 
pressure  of  21,000  kPa  (3000  psi) 

Dislance  from  center  of  source  to  center  of 
nearest  receiver  group 

50  m 

278  or  256  m 

Receiver  array 

48  groups  split  spread 

48  groups  bottom  drag  cable 

Group  interval 

100  m 

50  m  or  30.5  m 

Recording  method 

digital 

digital 

Sampling  rate 

2  ms 

2  ms 

Low  filters 

8  Hz 

5.3  Hz 

High  filters 

124  Hz 

128  Hz 

on  interpretation  error,  since  difficulties  in  seismic 
data  acquisition  and  analysis  may  occur  which 
cannot  be  evaluated.  For  example,  signal  identifi¬ 
cation  is  occasionally  a  problem.  Late-arriving  re¬ 
flections  can  be  misinterpreted  as  refraction 
events,  particularly  when  the  signal-to-noise  level 
is  low.  Another  problem  can  be  created  by  the 
need  to  assume  plane  layers  for  simplicity  in  inter¬ 
pretation,  whereas  localized  curvature  of  reflect¬ 
ing  and  refracting  surfaces  can  be  expected  due  to 
roughness  on  the  upper  surface  of  degrading  sub¬ 


sea  permafrost.  Both  of  these  problems  are  dealt 
with  primarily  by  looking  at  the  internal  consis¬ 
tency  of  the  results  where  there  is  duplication  in 
coverage,  and  finally  by  verification  with  indepen¬ 
dent  geophysical  techniques  and  drilling  results. 

In  addition  to  the  refraction  and  reflection  in¬ 
terpretations,  two  seismic  anomalies  were  mapped 
that  appear  to  be  directly  related  to  degrading  per¬ 
mafrost.  Anomalous  attenuation  rates  were  ob¬ 
served  and  natural  noise  sources  were  found  in 
some  areas.  The  occurrence  of  these  anomalies  in 
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Figure  2.  Sample  monitor  records  showing  high-velocity  refractors  interpreted  as  subsea  permafrost 


I 


v. 


Figure  3.  Ray  path  geometry  for  a  refractor  dipping 
at  an  angle  <t>.  The  shot  is  at  P  and  the  receiver  is  at  Q 
for  the  down-dip  shot.  The  shot  and  receiver  posi¬ 
tions  are  reversed  for  the  up-dip  shot. 


relation  to  distribution  of  permafrost  is  discussed 
in  the  text. 

The  process  of  attenuation  involves  the  absorp¬ 
tion  of  a  percentage  of  the  transmitted  energy  with 
each  cycle  of  the  wave  propagation.  Since  high 
frequencies  have  shorter  wavelengths  relative  to 
low  frequencies,  one  can  usually  perceive  attenua¬ 
tion  as  low  pass  filtering  or  pulse  broadening  of 
the  signal.  Small  attenuation  losses  are  usually  ob¬ 
served  in  saturated  marine  sediments  (Winkler  and 
Nur  1979).  However,  drastic  attenuation  can  be 
expected  when  a  liquid  and  gas  are  combined  in 
the  sediment  pores. 

Failure  of  materials  due  to  deformation  can 
cause  small-scale  seismic  events.  Some  movement 
and  differential  settlement  can  be  expected  in  the 
Harrison  Bay  sediments,  and  in  subsea  permafrost 
in  general,  as  permafrost  thaws  in  the  warmer  ma¬ 
rine  environment.  Local  induced  seismicity  has 
been  observed  in  other  regions  associated  with 
withdrawal  of  groundwater  and  extraction  of  fluid 
from  petroleum  reservoirs  (Yerkes  and  Castle 
1976).  Thaw  of  pore  ice  in  permafrost  can  also 
have  the  effect  of  reducing  the  pore  water  pres¬ 
sure,  creating  a  situation  analogous  to  fluid  ex¬ 
traction.  Localized  seismic  events  often  follow  dy¬ 
namite  blasts  in  permafrost,  most  likely  arising 
from  release  of  stress  developed  during  seasonal 
cooling  of  frozen  ice-rich  sediments  and  from  vol¬ 
ume  increases  due  to  ground  ice  formation.  Minor 
earthquake  activity  has  also  been  reported  in  the 
Beaufort  Sea  near  Barter  Island  (Barnes  and  Hop¬ 
kins  1978). 


METHODS 
Reading  records 

A  technique  for  reading  oil-industry  monitor 
records  has  been  developed  for  carrying  out  sub¬ 


sea  permafrost  surveys.  The  technique  was  origin¬ 
ally  developed  for  measuring  velocities  directly 
from  monitor  records  with  a  drafting  arm.  It  has 
now  been  adapted  for  use  with  a  digitizing  tablet 
for  a  small  computer. 

Three  types  of  waves  have  been  identified  on  the 
records  and  used  in  the  analysis:  refractions,  re¬ 
flections  and  surface  waves.  The  same  reading 
procedure  is  used  for  all  three  wave  types.  Each 
reading  consists  of  the  coordinates  (jc  and  t)  for 
the  tangent  point  on  the  time-distance  plot  plus  a 
slope  measurement  (c  =  dx/dt)  of  a  tangent  to  the 
curve.  This  information  is  con.erted  to  velocity 
data  and  depth  profiles  by  means  of  the  appropri¬ 
ate  equations  described  below. 

Refractions 

A  dipping  plane  layer  refraction  interpretation 
could  be  used  on  the  reversed  ice  shooting  records. 
Following  the  derivation  given  by  Grant  and  West 
(1965)  and  using  the  geometry  shown  in  Figure  3, 
the  critical  angle  is  given  by  te  =  sin'1  (F0/F,) 
where  V0  is  the  upper  layer  velocity  and  F,  is  the 
lower  layer  velocity.  The  apparent  velocity  in  the 
lower  layer  when  shooting  down-dip  is 

c'  =  Vq/ sin  (/e  +  <F)  =  K0/sin[sin"(F0/F,)  +  4>\ 

(1) 

where  <t>  is  the  dip  of  the  boundary.  For  shooting 
up-dip,  the  apparent  velocity  is 

c"  =  Ko/sin(/c-0)  =  Fg/sin  (sin-,(  F0/F,)  -  <£]. 

(2) 

Equations  1  and  2  can  be  combined  to  give  an  ex¬ 
pression  for  the  true  velocity  (F,)  in  the  lower 
layer: 

i/c'  +  l  /c*  =  sin(/'c  +  <t>)/  V0  +  sin(/c  -  <t>)/  V0 
=  2  sin/c  cos0/Fo 
=  2  cos <t>/Vv 
For  small  dip  angles, 

1/F,  =  '/2(l/r+  1/c*).  (3) 

Grant  and  West  show  that  the  distance  h,  from 
the  up-dip  shot  point  to  the  high-velocity  refrac¬ 
tor,  is 

ht  =(F0V2)[l-(F0/FJ)2],/i  (4) 


4 


where  /0  is  the  intercept  time  on  the  record.  We  did 
not  measure  the  intercept  time,  but  it  can  be  easily 
calculated  from  the  tangent  readings.  This  con¬ 
verts  eq  4  to 

h]  =(V0/2)(r-r/c-)[l-(V0/Vi)2r/2-  (5) 

For  purposes  of  constructing  seismic  cross  sec¬ 
tions,  the  depths  were  plotted  under  the  midpoints 
of  the  reversed  spread.  An  average  depth  ( h )  was 
calculated  for  an  array  based  on  values  from  its 
ends: 

h  =  ( F0/4)(r  +  /*  -  x-/c~  -  x'/e*) 

[i-(V^i  )2r,/2-  w 

Equations  3  and  6  are  the  required  equations  for 
making  velocity  and  depth  profiles  for  the  ice¬ 
shooting  data. 

There  are  no  reversed  profiles  for  the  marine 
survey  data;  therefore,  they  were  interpreted  as¬ 
suming  plane  horizontal  layers.  Equations  2  and  5 
can  be  rewritten  for  horizontal  layers  by  setting  <l> 
=  0  and  c  =  c*: 

V,  =  c  (7) 

h  =  ( V</2)(t -x/c)[\  -  ( V0/V{)\  (8) 

These  two  equations  allow  the  conversion  of  tan¬ 
gent  readings  from  the  marine  records  to  velocity 
and  depth  profiles  along  the  marine  lines. 

Reflections 

Reflection  data  analysis  was  based  on  assuming 
a  plane  horizontal  reflector  at  a  depth  h  under  a 
uniform  upper  layer  with  velocity  V0.  This  simple 
model  results  in  the  equation  of  a  hyperbola  for 
the  travel  time  t  and  the  distance  x  from  the  shot 
point  to  the  receiver  (see  Grant  and  West  1965): 

f£/2  =  j^  +  4/i2.  (9) 

Taking  differentials  on  both  sides  of  the  above 
equation  gives  an  expression  for  the  upper  layer 
velocity: 

IV^tdt  =  2j <dx.  (10) 

Since  dx/dl  =  c,  eq  10  can  be  arranged  to  give  the 
velocity  in  terms  of  the  tangent  readings: 

V0  =  (cx/t)v\ 


Combining  eq  9  and  1 1  gives  an  expression  for  the 
depth  to  the  reflector: 

h  =  (x/2)(ct/x  —  1 )  '/2 .  (12) 

Any  tangent  to  a  reflection  curve  can  be  con¬ 
verted  using  eq  1 1  and  12  into  a  depth  and  velocity 
determination  for  profile  construction  along  the 
shot  lines.  Where  possible,  a  number  of  tangents 
were  read  on  each  reflector,  so  that  scatter  on  the 
plotted  reflecting  horizons  could  help  to  indicate 
the  accuracy  of  the  determinations. 

Rayleigh  waves 

Based  on  the  surface  wave  or  Rayleigh  wave 
tangent  measurements,  profiles  were  made  of  the 
phase  velocity  along  the  survey  lines.  These  pro¬ 
files  proved  to  be  useful  for  identifying  the  signals 
from  the  small  quakes  or  noise  bursts  on  the 
records. 

Spatial  resolution 

The  spatial  resolution  of  the  data  from  this  type 
of  study  is  obviously  not  as  great  as  can  be  ob¬ 
tained  from  a  seismic  investigation  specifically  de¬ 
signed  to  study  offshore  permafrost.  For  refrac¬ 
tion  and  reflection  measurements,  a  number  of 
factors  must  be  considered,  including  geophone 
spacing,  signal  frequency,  and  complexity  of  the 
subsurface.  In  general,  the  horizontal  extent  of  a 
feature  that  can  be  detected  should  be  a  minimum 
of  three  detector  spacings.  This  means  that  the 
minimum  size  of  a  target  that  can  be  resolved  is 
around  300  m  for  the  ice-shooting  data,  and 
around  150  m  for  the  marine  survey  data.  The 
minimum  vertical  thickness  of  a  detectable  high- 
velocity  layer  is  determined  by  the  wavelength  of 
the  refracted  signal  (Sherwood  1967).  Resolution 
is  possible  to  approximately  Vi  wavelength  or 
about  50  m  for  these  data.  In  addition,  Sher¬ 
wood’s  results  show  that  thin  layers  (less  than  30 
m  thick)  might  be  observed  at  shallow  depths,  but 
the  signals  from  these  would  be  in  the  form  of 
plate  waves  at  a  reduced  velocity  and  amplitude. 

A  simplifying  assumption  was  used  for  the  re¬ 
fraction  depth  determinations.  The  upper  layer 
velocity  was  taken  as  1.8  km/s  for  all  profiles. 
This  means  that  the  water  layer,  0  to  16  m  deep, 
was  combined  with  the  low-velocity  bottom  sedi¬ 
ments  to  make  a  single  upper  layer.  Upper-layer 
velocities  were  observed  to  range  from  1.6  to  2.0 
km/s.  Therefore,  the  error  introduced  by  assum¬ 
ing  1 .8  km/s  could  be  as  much  as  30%  under  rare 
circumstances  (see  Appendix  A). 
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Figure  4.  Velocity  data  from  refraction  analysis  of 
Harrison  Bay  records. 


Refraction  velocities  and  depth  determir  .dons 
from  single-ended  marine  records  are  subject  to 
errors  caused  by  dipping  layers.  Our  interpreta¬ 
tions  indicate  that  dips  are  normally  less  than  3%. 
The  corresponding  maximum  error  is  approxi¬ 
mately  5%  in  velocity  measurements  and  2%  in 
depth  determinations  (see  Appendix  A). 

The  assumption  of  horizontal  layers  for  the  re¬ 
flection  interpretation  does  not  result  in  signifi¬ 
cant  errors.  Error  calculations  in  Appendix  A 
show  that  a  3%  dip  usually  results  in  a  1%  error  in 
velocity  and  2%  in  depth. 

Anomalies 

Examples  of  the  anomalous  features  observed 
on  the  records  are  shown  in  the  text,  with  their  dis¬ 
tribution  indicated  by  anomaly  maps.  The  attenu¬ 
ation  map  was  based  on  a  simple  qualitative  pro¬ 
cedure.  The  records  examined  normally  have  a 
peak  frequency  of  30  to  40  Hz.  All  records  with 
peak  frequencies  lower  than  1 5  Hz  were  classed  as 
having  been  subject  to  considerable  natural  filter¬ 
ing  and  were  grouped  in  the  attenuation  zones. 

No  formal  attempt  was  made  to  determine  the 
attenuation  coefficients  for  comparison  with  exist¬ 
ing  data  on  marine  sediments.  Such  a  procedure 
would  require  spectral  analysis  to  establish  relative 
amplitudes  at  the  different  frequenc  es.  This  was 
not  possible  with  the  available  records  and  equip¬ 
ment. 


Examples  of  natural  seismicity  are  presented, 
and  an  anomaly  map  indicating  the  distribution 
was  prepared  for  comparison  with  other  maps  in¬ 
cluding  the  distribution  of  high-velocity  material 
and  attenuation  anomaly  distribution.  Compara¬ 
tive  profiles  were  constructed  of  phase  velocities 
from  the  natural  seismicity  and  phase  velocities 
from  surface  waves  generated  by  the  air  guns.  This 
velocity  comparison  was  helpful  for  wave  iden¬ 
tification. 


RESULTS  AND  DISCUSSION 

Seismic  velocity  distribution 

The  range  of  seismic  velocities  observed  in  the 
Harrison  Bay  region  is  shown  in  Figure  4.  The  low 
velocity  peak  is  in  the  range  usually  associated 
with  thawed  material,  and  the  higher  velocities 
(>2.0  km/s)  are  usually  considered  representative 
of  chemically  bonded  or  ice-indurated  material. 

The  cross  sections  in  Figures  5  and  6  represent 
the  two  velocity  regimes  found  in  Harrison  Bay. 
Both  lines  have  a  high-velocity  structure  extending 
approximately  25  km  offshore.  The  eastern  line 
(Fig.  5)  shows  a  systematic  thickening  of  the  first 
layer  (low-velocity  layer)  with  increasing  distance 
from  shore.  The  mean  velocity  of  this  upper  layer 
(1.8  km/s)  falls  in  the  range  that  represents  little 
or  no  ice-bonding  of  the  sediments.  The  second- 
layer  velocities  of  3  to  4  km/s  are  consistent  with 
ice-bonded  material.  In  addition,  the  deeper  laye- 
shows  a  velocity  decrease  with  depth  (Fig.  5b), 
suggesting  that  the  materials  become  warmer  with 
depth.  This  thickening  of  the  low-velocity  layer 
with  increasing  distance  from  shore  and  the  de¬ 
crease  in  velocity  of  the  second  layer  with  depth 
support  a  model  of  degrading  ice-bonded  perma¬ 
frost  in  an  area  of  active  marine  transgression. 
The  depth  of  the  lower  layer  is  greater  than  the 
depth  of  a  similar  layer  observed  near  Prudhoe 
Bay,  which  indicates  differing  geological  and  per¬ 
mafrost  conditions  in  the  two  regions. 

The  second  set  of  profiles  (Fig.  6)  illustrates  the 
greater  complexity  encountered  in  the  western  half 
of  Harrison  Bay.  Four  separate  velocity  zones  are 
encountered  along  this  line.  The  first  segment,  in¬ 
cluding  the  onshore  records  and  several  offshore 
records,  is  anomalous  because  of  the  low  average 
velocities  between  the  surface  and  the  deep  reflec¬ 
tors.  A  second  segment  near  the  shore  has  a  shal¬ 
low,  high-velocity  refractor.  The  remaining  off¬ 
shore  half  of  the  line  has  two  distinct  zones,  one 
deep  and  one  shallow,  that  may  be  partially  bond¬ 
ed  or  may  represent  a  change  in  material  type. 

The  deep  reflections  at  the  south  end  of  the  line 
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Figure  5.  Profiles  from  the  eastern  part  of  Harrison  Bay. 
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Figure  6.  Profile  from  the  western  part  of  Harrison  Bay. 


(Fig.  6),  and  also  on  line  D32  (see  Appendix  C), 
are  surprising  because  there  is  only  a  thin  layer  of 
high-velocity  material  at  the  surface.  Under  that 
lies  a  thick  section  with  low  average  velocity.  The 
low-velocity  zone  could  be  explained  as  a  thaw 
zone  formed  beneath  a  lake  that  historically  occu¬ 
pied  the  surface  in  this  area.  This  would  imply  that 


recent  near-surface  temperatures  are  low  due  to 
exposure  of  the  surface  as  a  result  of  draining  of 
the  thaw  lake  basin.  The  deeper  sediments  down 
to  400  m  are  warmer,  reflecting  the  previous 
warming  cycle  induced  by  the  lake  environment. 
The  low  velocities  could  also  be  explained  by  a 
contrast  in  material  properties.  If  the  thin  surface 
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Figure  7.  Distribution  of  high-velocity  material  in  Harrison  Bay.  A  continuous  high  velocity  is  present  be¬ 
tween  the  shore  and  the  first  contour  line.  Reflection  readings  show  that  this  feature  can  be  followed  farther 
offshore  to  the  second  contour  line  before  it  disappears. 


Figure  8.  Depth  to  the  top  of  the  high  velocity  refractor. 


layer  were  ice-bonded  coarse-grained  material, 
and  the  rest  of  the  section  were  fine-grained  with 
less  complete  ice-bonding,  the  usual  thick  layer  of 
permafrost  (300-600  m)  could  be  present. 

Velocity  profiles  and  cross  sections  like  those  in 
Figures  3  and  6  were  constructed  for  all  the  shot 


lines  in  Harrison  Bay  (Appendices  B  and  C). 
These  data  suggest  the  distribution  of  ice-bonded 
permafrost  when  compiled  on  a  set  of  area  maps. 
The  structural  zones  are  shown  in  Figure  7.  A  two- 
layer  situation  exists  in  the  zone  nearest  shore, 
with  a  high-velocity  layer  at  depth  that  decreases 


Figure  9.  The  refraction  velocity  contours  for  the  high  velocity  layer. 


in  velocity  and  increases  in  depth  with  distance 
from  shore.  Beyond  this  zone,  the  continuation  of 
the  high-velocity  structure  is  suggested  by  a  deep 
reflector.  Further  offshore,  there  are  no  velocities 
above  2.0  km/s,  making  it  difficult  to  determine 
whether  the  materials  are  ice-bonded.  If  they  are, 
the  low  velocity  would  indicate  that  they  are  warm 
or  have  relatively  low  ice  content. 

A  contour  map  of  the  depth  to  the  refractor 
(Fig.  8)  reveals  that  the  high-velocity  layer  is  at  or 
near  the  surface  at  the  shore  and  dips  to  the  north¬ 
east  under  the  bay.  The  steepest  dips,  approxi¬ 
mately  3%,  are  near  shore  off  Atiguru  Point  and 
just  west  of  the  Colville  Delta.  Lesser  slopes,  ap¬ 
proximately  1%,  are  seen  further  offshore  and 
near  the  Cape  Halkett  coastline. 

The  shallow  depth  of  the  refractor  observed  in 
the  western  part  of  the  bay  corresponds  with  ob¬ 
servations  made  to  the  west  near  Lonely,  which  is 
in  the  same  geological  setting.  Harrison  and  Oster- 
kamp  (1981)  established  a  probe  line  offshore 
from  Lonely  and  observed  shallow  ice-bonded 
permafrost  8  to  15  m  below  the  seabed  out  to  at 
least  7.8  km  from  shore. 

A  contour  map  of  refractor  velocities  (Fig.  9) 
also  helps  to  illustrate  the  difference  between  the 
eastern  and  western  parts  of  Harrison  Bay.  Veloci¬ 
ties  west  of  Atiguru  Point  are  more  variable  and 
decrease  more  rapidly  with  distance  from  shore 
than  velocities  east  of  the  point.  These  differences 
suggest  that  a  contrast  must  exist  between  the  ma¬ 
terial  types  and  also  between  the  geological  histo¬ 


ries  of  the  two  regions.  The  western  part  of  the 
bay  may  have  been  an  extension  of  the  low,  lake- 
covered  coastal  plain  found  west  of  Cape  Halkett. 
The  shallow  depth  of  the  refractor  in  this  area  and 
the  noticeable  decrease  in  velocity  with  distance 
from  shore  would  be  expected  in  an  area  of  fine¬ 
grained  material  subject  to  active  coastal  retreat. 

Attenuation 

Figure  10a  is  a  record  which  illustrates  the  sig¬ 
nature  of  the  air  guns  and  recorder  system.  It 
shows  an  attenuation  rate  in  the  sediments  which 
we  do  not  consider  anomalous.  The  dominant  fre¬ 
quency  of  the  traces  0.28  km  from  the  source  is  in 
the  20-  to  40-Hz  range.  At  the  far  end  of  the  array 
(2.6  km  from  the  source)  the  dominant  frequency 
is  approximately  20  Hz.  In  addition,  strong  high- 
frequency  modes  of  80  Hz  can  be  seen  as  later  arri¬ 
vals  on  these  far  traces.  Figure  10b  shows  an  ano¬ 
malous  situation  where  high  frequencies  are  pres¬ 
ent  near  the  source.  Most  noticeable  is  the  strong 
mode  at  about  80  Hz.  The  dominant  signal  at  the 
far  end  of  this  array  is  10  Hz.  A  more  common 
type  of  anomalous  record  is  illustrated  in  Figure 
10c.  The  entire  record,  excluding  background 
noise,  is  composed  of  frequencies  of  less  than  15 
Hz.  Similar  attenuation  is  seen  on  the  ice-shooting 
records  (Fig.  lOd  and  lOe).  We  attribute  these 
anomalies  in  both  the  ice  and  marine  records  to  at¬ 
tenuation  in  the  sediments  for  a  number  of  rea¬ 
sons.  First,  the  field  notes  do  not  indicate  any  sig¬ 
nificant  air  gun  malfunctions,  and  no  changes  in 
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Figure  10.  Sample  records  for  the  attenuation  stuay. 


cord  which  has  a  dominant  frequency  of  less  than  IS 
whole  spread. 
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Figure  //.  Distribution  map  of  the  " shallow  gas”  anomaly.  The  high  frequency  part  of  the  seismic  signals 
has  been  naturally  filtered  in  the  shaded  zones. 


the  recording  filters  were  made.  Secondly,  these 
anomalies  can  be  traced  from  line  to  line  within  a 
survey.  Finally,  the  near-shore  extent  of  the  atten¬ 
uation  zones  defined  by  both  the  marine  and  ice 
shooting  surveys  are  in  agreement.  This  seems  to 
be  fairly  positive  evidence  that  the  anomalies  are 
related  to  seabed  characteristics,  since  the  surveys 
had  different  energy  sources  and  were  conducted 
during  different  times  of  the  year. 

The  attenuation  effect  does  not  appear  to  di¬ 
minish  the  amplitude  of  the  low-frequency  end  of 
the  spectrum  for  records  such  as  Figure  10c.  This 
can  be  explained  if  these  signals  constitute  the  low¬ 
est  mode  in  the  near-surface  waveguide  which 
passes  the  recorder  filter  range.  There  would  be  no 
significant  energy  recorded  at  any  lower  frequen¬ 
cies,  and  these  signals  would  be  maintained  at 
nearly  constant  amplitudes  across  the  record  by 
the  automatic  gain  control  of  the  receiver  ampli¬ 
fiers. 

Figure  1 1  outlines  the  areas  where  the  dominant 
frequency  of  the  reflected  or  refracted  signals  was 
reduced  from  approximately  30  Hz  to  less  than  IS 
Hz.  These  areas  are  of  considerable  interest  be¬ 
cause  strong  attenuation  has  been  associated  with 
the  presence  of  gas  in  the  pores  of  sediments.  Lab¬ 
oratory  studies  by  Winkler  and  Nur  (1979)  have 
shown  that  pore  gas  concentrations  as  low  as  a  few 
percent  by  volume  can  cause  attenuation  to  in¬ 
crease  by  a  factor  of  5  over  water-saturated  val¬ 


ues.  Therefore,  the  high  attenuation  observed  on 
Harrison  Bay  records  could  be  explained  by  the 
presence  of  natural  gas  in  the  sediments.  Shallow 
gas  has  been  reported  near  Prudhoe  Bay  by  Bou¬ 
cher  et  al.  (1981),  based  on  the  results  of  high  reso¬ 
lution  seismic  studies.  It  has  also  been  observed  in 
the  Canadian  Beaufort  Sea  off  the  Mackenzie  Del¬ 
ta  by  Neave  et  al.  (1978),  based  on  evidence  from 
drilling  and  seismic  studies.  As  a  result,  we  feel 
that  there  is  considerable  evidence  that  the  attenu¬ 
ation  zones  in  Figure  1 1  can  be  viewed  as  the  shal¬ 
low  gas  distribution  for  Harrison  Bay. 

The  depth  of  the  gas  deposits  has  not  been  de¬ 
termined.  However,  the  attenuated  signals  from 
both  refraction  and  reflection  analysis  penetrate 
to  depths  in  the  range  20  to  400  m.  Consequently, 
the  gas  deposits  must  lie  in  this  range. 

According  to  the  stability  conditions  for  gas  hy¬ 
drates,  a  solid  form  of  natural  gas  and  water  can 
be  found  in  permafrost  at  a  depth  of  100  m  or 
greater,  depending  on  the  composition  of  the  gas 
(Davidson  et  al.  1978,  APOA  1978).  Therefore, 
there  could  be  some  gas  in  hydrate  form  within 
and  below  the  ice-bonded  layer  in  these  areas  of 
shallow  gas. 

The  proposed  gassy  zone  in  Figure  1 1  covers  a 
large  portion  of  the  area  in  Figure  9,  where  ap¬ 
parently  active  degradation  of  ice-bonded  perma¬ 
frost  is  taking  place,  suggesting  a  direct  relation¬ 
ship  between  the  two  zones.  The  source  of  the  gas 
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Figure  I 
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Figure  13.  The  distribution  of  natural  seismic  activity  in  Harrison  Bay.  The  strongest  events  and  highest  ac¬ 
tivity  are  in  a  band  close  to  shore. 


is  unknown.  It  may  be  biogenic  gas  formed  in 
place  or  it  may  be  from  depth,  freed  by  thawing  of 
impermeable  zones,  or  from  decomposition  of  hy¬ 
drates. 

Low-level  natural  seismicity 

Most  of  the  noise  on  the  records  examined  as 
part  of  this  study  seems  to  be  a  form  of  natural 
seismicity  associated  with  areas  of  shallow  subsea 
permafrost.  The  best  examples  were  observed  in 
the  Sagavanirktok  Delta  and  Tigvariak  Island  area 
just  east  of  Prudhoe  Bay.  Similar  events  were 
found  in  Harrison  Bay.  Figures  I2a-d  reproduce 
four  of  the  clearest  events  from  the  records  east  of 
Prudhoe.  They  are  not  found  at  any  particular 
place  or  time  on  the  records.  The  Harrison  Bay 
events  (Fig.  12e-g)  are  similar,  aside  from  the  fact 
that  they  have  more  dispersion  and  low  ampli¬ 
tudes.  A  good  correlation  of  the  activity  turns  out 
to  be  with  the  distribution  of  ice-bonded  perma¬ 
frost.  However,  the  prospect  that  this  noise  may 
be  generated  by  some  other  feature,  such  as  ice 
movement,  cannot  be  discounted. 

Regions  in  which  seismic  events  were  observed 
are  shown  in  Figure  13.  These  regions  correspond 
with  both  the  location  of  the  high-velocity  mater¬ 
ials  and  the  regions  of  shallow  gas.  The  noise 
could  be  related  to  energy  released  during  perma¬ 
frost  degradation.  In  such  a  process,  the  energy 
could  originate  from  expanding  gas  or  adjustment 


of  seabed  material  to  differential  settlement 
caused  by  permafrost  thaw. 

Figure  14  compares  the  phase  velocity  of  Ray¬ 
leigh  waves  from  the  air  guns  with  the  phase  velo¬ 
cities  of  the  induced  seismicity  for  Line  311.  They 
both  lie  in  the  same  velocity  range  (0.3  to  1.1 
km/s),  which  is  below  the  compressional  wave  ve¬ 
locity  in  water  (1.4  km/s)  and  also  much  less  than 
the  plate  wave  velocity  in  sea  ice  (2.5-3.0  km/s). 
This  observation  helps  to  confirm  that  the  noise  is 
associated  with  the  sediments  and  not  with  waves 
or  ice.  The  group  velocities  (not  shown)  are  less 
than  the  phase  velocities,  so  their  dispersion  is  like 
Rayleigh  wave  normal  modes.  Therefore,  the 
measured  phase  velocities  lie  between  the  Rayleigh 
wave  velocities  of  the  thawed  layer  and  the  frozen 
layer.  See  Ewing  et  al.  (1957),  for  a  discussion  of 
these  velocities. 

The  amplitudes  of  these  signals  cannot  be  deter¬ 
mined  due  to  automatic  gain  control  on  the  re¬ 
cording  equipment.  However,  they  appear  to  have 
an  energy  output  similar  to  that  of  the  air  guns. 
Because  there  are  no  recordings  without  air-gun 
shots,  it  is  not  clear  whether  or  not  the  noise  was 
triggered  by  the  shots,  like  the  aftershocks  that  oc¬ 
cur  in  land  seismic  data  in  permafrost.  In  any 
case,  the  noise  suggests  that  unstable  conditions 
are  present  in  areas  of  actively  degrading  perma¬ 
frost. 
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Figure  14.  Phase  velocities  from  natural  signals  compared  to  the 
phase  velocity  of  surface  waves  from  the  air  guns  for  Line  311. 


SUMMARY 

Velocity  profiles  for  Harrison  Bay  show  that 
high-velocity  material  interpreted  as  ice-bonded 
permafrost  is  common  in  this  region.  The  high- 
velocity  ice-bonded  sediments  extend  at  least  25 
km  from  shore.  There  is  a  considerable  variation 
in  the  apparent  permafrost  distribution,  which  in¬ 
dicates  the  significant  effect  of  material  types, 
past  inundation  and  geological  history  on  the  rate 
of  subsea  permafrost  degradation.  In  particular, 
the  part  of  the  bay  west  of  Atiguru  Point  is  more 
chaotic  in  terms  of  velocity  and  depth  variations. 
The  eastern  half  of  the  bay  has  far  more  regular 
permafrost  profiles.  Across  the  whole  northern 
edge  of  the  study  area,  it  appears  that  thermal  de¬ 
gradation  has  modified  the  bonded  permafrost  to 
depths  as  great  as  300  m. 

The  anomalous  attenuation  on  the  seismic  rec¬ 
ords  observed  in  Harrison  Bay  was  interpreted  as 
an  indicator  of  free  gas  in  the  section  above  the 
frozen  sediments.  This  means  that  gas  hydrates 
can  occur  at  greater  depths  and  may  be  the  source 
of  the  free  gas  liberated  as  permafrost  thawed. 

The  swarms  of  small  quakes  or  seismic  noise  ob¬ 
served  in  the  nearshore  areas  were  interpreted  as 
energy  released  due  to  thermal  modification  of 
this  recently  inundated  permafrost. 
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APPENDIX  A:  ERROR  ESTIMATES 

The  first  error  estimate  arises  from  using  an  average  upper  layer  velocity  of  1 .8  km/s  for  refrac¬ 
tion  calculations  instead  of  the  local  value.  A  worst  case  situation  is  examined  to  illustrate  how 
much  error  is  introduced.  The  remaining  error  problems  involve  dipping  layers  when  the  interpre¬ 
tation  assumes  horizontal  layers.  For  these  cases,  we  used  typical  readings  and  calculated  the  dif¬ 
ference  between  the  velocity  and  depth  estimates  for  horizontal  layers  and  layers  with  a  3%  dip. 

When  the  correct  local  velocity  K,  is  used  in  eq  8  of  the  Methods  section,  the  depth  is: 

ht=(t~xlc)(VJ2)  [1  -(Kl/K,)2]-*  .  (13) 

When  the  average  velocity  VA  ~  1 .8  km/s  is  used  instead,  the  erroneous  depth  estimate  is 


hA=(t~xlc)(VA/2 )  [1  -(VJV,)*)-*  . 


The  proportional  error  in  the  depth  estimate  is 


(*L -*a)/Al  =  1 -*a/*l 


^a  n-(*w2i* 

A  Li-OV^i)2] 


This  expression  gives  the  largest  error  when  the  local  velocity  is  large  and  the  lower  layer  velocity 
Vt  is  small.  From  the  velocity  histogram  in  Figure  4,  we  can  choose  VL  -  2.0  km/s  and  K,  =  2.3 
km/s.  This  highly  unlikely  combination  of  velocities  results  in  an  error  of  30%. 

The  effects  of  dip  on  the  interpretation  of  apparent  lower  layer  velocity  from  single-ended  re¬ 
fraction  data  can  be  found  from  eq  1  : 

c=F0/sin[sin-«(ro/F1)+^J  .  (16) 

The  real  veloicty  in  the  lower  layer  can  be  found  by  solving  this  equation  for  Vl : 

Vx  =  VQ/sin  [sin-1  (VQ/c)-<t>]  .  (17) 

The  proportional  error  in  the  lower  velocity  from  using  a  horizontal  layer  model  is 

(V i  ~c)/Vl  =l-c/Vl 


=  1  -(c/V0)  sin  [sin-1  ( V0/c)-<p J  . 


The  estimated  depth  for  the  horizontal  layer  model  is  given  by  eq  8: 


hH=(V0/2)(t-x/c)  [l-(F0/c)*j-* 


This  equation  can  have  K,  replaced  by  a  substitution  from  eq  17.  Then  eq  19  and  20  can  be  used 
to  find  the  proportional  error  in  the  depth: 


(*»d  "*h)^d  s  1  -  cos  [sin-1  (V0lc)-4>\  fl  -(K0/c)J ]'*  . 


Using  a  typical  set  of  velocities,  V0  =  1 .8  km/s  and  c  =  3.66  km/s.  A  slope  of  3%  results  in  a  2% 
error  in  depth  according  to  eq  21 ,  and  the  corresponding  error  in  the  velocity  is  5%  from  eq  18. 

For  the  reflections  from  a  plane  boundary  which  dips  at  an  angle  0,  the  arrivals  on  the  record 
still  form  a  hyperbolic  curve;  however,  the  hyperbola  is  not  centered  with  respect  to  the  shot  point. 
The  arrival  times  can  be  calculated  from  an  image  source  R  (Fig.  Al)  which  is  at  a  depth  ?  =  2 h 
cos 0  and  displaced  by  a  horizontal  distance  £  =  2h  sin0  from  the  true  source  P. 


Figure  Al.  Ray  path  geometry  for  a  reflector 
dipping  at  an  angle  0. 


The  travel  time  equation  is 
f2  -  (*  +  £)2  +  cos 0)2 

n 

t 2  =  (1/F2)  (4h2  cos2  0  +  4 h2  sin2  0  +  4 hx  sin0) 

t 2  V2  =4/r2  +x2  +4Axsin0.  (21) 

Taking  differentials  on  both  sides  of  the  equation  gives 
2 1  V2dt  =  2x dx  +  4h  sin0  dx 

V2  =  (x/t)  ( dxldt )  +  (2h  sin 0/r)  (dx/dt) .  (22) 

The  estimated  velocity  for  the  horizontal  layer  interpretation  was 


This  can  be  substituted  into  eq  23: 

V2  =  F*,  +  2h  sin 0  vyx 

V0  =  VH  [1  +  2/r  Sin0/jc]  ^  .  (24) 


The  proportional  error  in  velocity  is 

(V0-VH)IV0  =  \-[\  +  2A  Sin0/jcl . 
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(25) 


Now  compare  the  horizontal  layer  depth  estimate  in  eq  12 
hH=(xll)  [(ct/x)-l]* 
to  the  distance  to  reflector  in  eq  22 
t2  V2  =  Ah2  +  x2  +  Ah  x  sin#  . 

To  eliminate  we  use  eq  24  and  25: 

ctx  [1  +  (2 h  sin#/x)]  =  Ah2  +  x2  +  4 hx  sin#  . 

This  is  rearranged  as  a  quadratic  equation  with  h  as  the  unknown: 

Ah2  +  (4x  sin#  -  2ct  sin#)  h  +  x2  -  ctx  =  0  . 

Solving  for  h  gives 

hd  =  (1/8)  (2ct  - 4x)  sin#  ±  (1/8)  £(2ct  - 4x)2  sin2#  -  16  (x2  -  ctx)J* 

=  [(ct/4)  -  (x/2)]  sin#  ±  (1/4)  [(ct  -  2x)2  sin2#  -  4(x2  -  ctx)] 'A 

The  relative  error  in  depth  is 

hd~hH  _  J _ (X/2)  [ (ct/x)  -  1  ]  Vi _ 

hd  1  /4  (ct  -  2x)  sin#  ±  1  /4  J(ct  -  2x)2  sin2#  -  4x  (x  -  ct)J  Vl 

For  a  typical  reflection  reading,  c-  2.0  km/s,  x  -  l  .0  km  and  t  =  0.63  s.  With  a  dip  of  3%  on 
the  reflector,  we  find  a  1%  error  in  the  velocity  determination  and  a  2%  error  in  the  depth. 


(26) 


(27) 


(28) 


21 


APPENDIX  B:  VELOCITY  PROFILES 


The  following  profiles  are  plotted  with  south  on  the  left.  The  velocity  decreases  consistently 
from  left  to  right  except  on  Line  303,  which  is  parallel  to  the  strike  of  the  structure.  Where  there 
are  two  profiles  the  low-velocity  profile  is  the  near-surface  velocity.  In  two  cases,  Lines  307  and 
31 1 ,  the  velocity  rises  at  the  right  end  of  the  line.  The  right  ends  of  these  lines  approach  the  coast 
near  Cape  Halkett. 

Both  refraction  and  reflection  velocities  are  plotted  on  most  of  the  following  lines.  Refraction 
data  are  indicated  by  (+)  and  reflection  data  by  (0). 

The  profiles  are  identified  by  a  line  number.  The  location  of  these  lines  in  Harrison  Bay  can  be 
found  in  Figure  1 . 
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APPENDIX  C:  SEISMIC  CROSS  SECTIONS 


The  following  depth  profiles  show  a  northeast  dip  of  1-3%  for  the  refracting  horizon.  Three 
lines  (D26,  D30  and  D31)  have  the  high  velocity  material  rising  to  meet  the  shoreline.  Refraction 
multiples  are  found  on  a  number  of  profiles.  Vertical  expression  is  28: 1 . 

The  following  legend  applies  to  the  profiles: 

+  Refraction  depth  (velocity  <  2.0  km/s) 

*  Refraction  depth  (velocity  >  2.0  km/s) 

0  Reflection  depth. 

The  location  of  the  lines  can  be  found  in  Figure  1 . 
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